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ROBUSTIFYING THE KALMAN FILTER; 

Protection Against Symmetrically Straggling Measurement Errors 



D. P. Gaver 
P. A. Jacobs 



1. INTRODUCTION. 

Tracking and prediction algorithms based on simple Gaussian (normal 
distribution) measurement errors and structural models are commonly used in 
practice under the name of KALMAN Filters. If (a) measurement errors are not 
suitably Gaussian, e.g., if occasional outliers occur or (b) true structural 
behavior is not simple, perhaps displaying apparently discontinuous behavior 
caused by unfavorable sensor-target orientation, then traditional filter 
performance may dramatically degrade. In this paper, we will propose and 
study procedures based on an elaborated model of die KALMAN-type but with 
the measurement errors coming from a family of possibly suitable non-Gaussian 
distributions (e.g., Student-t) to represent, and suitably compensate for more- 
thick-tailed-than-Gaussian measurement error, i.e., distributions with long 
straggling tails having the tendency to produce symmetric outliers. 

In particular the basic stochastic model considered here is 



6n - 6n-l + 0>n 


(i.i) 


Yn = 9n + En 


(1.2) 



where { co n } are independent normal/Gaussian random variables with mean 0 
and variances {T n } and { e n ) are independent random variables having mean 0. 
The random variable 0 n is unobservable. The random variable Y n is interpreted 
as the observation of 0 n made with measurement error e n ; £n is not Gaussian, 

but controllably long-tailed. The problem is to estimate 0 n from Yj,..., Y n in 
the simple recursive fashion that characterizes the classical KALMAN filter. 
Expression (1.1) is a simple random walk and does not represent very 
interesting dynamics, but does provide suggestive illustrations. 

In the next, or second, section, we will describe a procedure, the ALMA 
(standing for KALMAN with outliers suppressed), which is based on a model in 
which the components of the error sequence (e n ) have a Student-t distribution. 



In the third section, the traditional KALMAN procedure will be 
described. It is based on the assumption that components of { e n } have iid 
normal distributions. Finally, a robust procedure due to West [1981 1 will be 
described. 

In section 4 results of an extensive simulation experiment will be 
presented and discussed. The simulation experiment compares the various 
procedures. The results indicate that the ALMA procedure is significantly 
better than the KALMAN when the true measurement error distribution is 
Student-t. Further, there is not much lost in using the ALMA procedure instead 
of the KALMAN when the true measurement error distribution is normal. 

2. THE ALMA FILTER AND RELATED PROCEDURES. 

While many measurement errors of physical quantities are approximately 
normal, especially "in the middle" of their distribution, there can well be 
thicker-than-normal/Gauss tails and also occasional extreme outliers; that these 
can have seriously degrading effects in regression-like problems has been the 
subject of considerable research; we cite books by Mosteller and Tukey (1977), 
Huber (1981), Hampel (1986); in the time-series context the article by Martin and 
Yohai (1986), which contains many references; also lately the articles by West 
and his associates (1981,1985); it is to West's approach that our methodology 
should best be compared. 

One way to model these features is to extend the tails of the normal by 
continuous scale mixing. Such an approach can lead to the Student-t form, and 

to many other useful forms as well. We will assume here that { e n } are 
independent random variables, now having in the Student-t distribution with 
mean 0, scale a n ( not the standard deviation) and d degrees of freedom; that is, 




( 2 . 1 ) 



Let yj denote the i* measurement and y n = (yj, ..., y n ). Assume that 9 n - 
ll y n_ l has a normal distribution with mean m n .j and variance C n -i. Since co n is 
assumed to have a normal distribution with variance x n , 0 n l y 0 * 1 has a normal 



distribution with mean m n -i and variance Cn=C n -i+T n . Thus, from (1.1), (1.2), 
and (2.1) 



P {0 n ed0,Y n €dy |Yi=y 1 ,...,Y n . 1 =y n . 1 } 




n 



( 2 . 2 ) 




where the approximation replaces the expression in the exponent by an 
approximating quadratic in 0. 

2.1 The ALMA Procedure. 

The ALMA procedure provides a Gaussian approximation to the 
distribution of 0 n | y n , but one that emphatically differs from the classical linear- 
in-observations form. Following an argument in Gaver et al. [1986], 
differentiate both sides of (2.2) with respect to 0 to obtain 



e-p(y) e-m,- 1 [ d+1 0-y 



C(y) r # d 2 

n 




(2.3) 



Equating the terms involving 0 results in the following equation: 



where the weight 



, , d+1 

w(y) = — 



(2.5) 



1 + 



<0-y\ 



n y 



1 

a 



Furthermore, equating the constant terms results in 



F(y) m n .i 



(y) c* 



+ w 



(y>4- 



(2.6) 



r rhe ALMA procedure approximates 0 n I y n by the normal distribution 
having mean 



m n = S, 

and variance 



m n- 1 h + 

u n 



(2.7) 



C n = 



1 , , 1 
~T~ + w(y n ) -r 

:: ot 



(2.8) 



where 



w (y n )= 



1 + 



fe-yrA 



°n / 



(2.9) 



1_ 

d 



Note that the weight w(y n ) involves the unknown 0. One implementation 
uses approximate weights of the form 



w k(y n ) = 



l 



('n~ m n-l > \ k 



(2.10) 



1 + 



When k=l, m n _i is used in place of 0 in (2.9). 
When k=j, 0.5(m n .i +y n ) is used in place of 0. 



The basic ALMA procedure is to evaluate wj^yn) and then use it to find 



C„ = 



' 1 . , 1 1 
7# + w k(y n ) ~2 

n a n 



-1 



and 



m n — ^-n 



^ m n-l w k (y n )y n ^ 
C* o 2 



( 2 . 11 ) 



( 2 . 12 ) 



The point estimate of 0 n given y n is 0 n =m n and an estimate of the variance 

of 0 n is C n . Thus the procedure provides a particular Gaussian posterior 
approximation. In other similar contexts, non-linear filters for example, it has 
been suggested that the procedure (2.10) - (2.12) be iterated with the newly- 
computed m n , replacing m n _i in (2.10) - (2.12) in each iteration. In the 
simulations 0, 1 and 2 iterations were implemented, and the results compared. 

2.2 The Biweight. 

The ALMA procedure is an iterative reweighting procedure, hi the 
ordinary regression context another weight has been suggested: the so-called 
(Tukey) biweight, cf. Mosteller and Tukey (1977). In our context, the biweight 
procedure can replace the weight w^(y) in the ALMA procedure with the 
biweight 



wfi(y)= 







f 


'*? 


-|2 


f 


-1\2 


< 


1 - k 


(y-m n -J ac n \/tt 




if k 


(y-mn-ii aa n A/ry 








V L J 


/ 


a 




/ 


< 




0 






otherwise. 





<1 



(2.13) 



The variance of a Student-t distribution with d degrees of freedom and 

2 d 

scale a is crg^if d>3, otherwise being infinite. Hence the (bi)weight w B (y) uses 
the measurement y if I y I is within a standard deviations of m n _i, the estimate of 
0 n .j. The weight is zero if the deviation is greater. 

As was done in the basic ALMA procedure, 0, 1, and 2 iterations of 
(2.10)-{2.12) were tried, with w B (y n ) replacing Wk(y n ), for values of a=5, 7, 9 

and k=l , 0.25. 



2.3 Aspects of the Likelihood Procedure . 

It is possible for the likelihood function (2.2) to exhibit two local 0- 
maxima. In such a case, the likelihood procedure approximates the local 
maxima and chooses the one which globally maximizes the likelihood. 

To examine the details let 
f(0) = — -InP { 0 e d0,Y edy) 

HA n n 3 



Now it is clearly possible for f(0)=O to have multiple roots. To be 
specific, f(0)=O for those 0 satisfying 



The properties of this cubic-in-0 equation can be deduced from classical 
results. 




(2.14) 



0 = 0 3 +0 2 (-2y-m ,) 

J n-1 

+ 0^d + y 2 +(d+l)C* + 2ym nl j 

+ r- m n-l^ d - m n.iy 2 ' (d+1) y C nT 



(2.15) 



Let 
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r 

C # , i 


fc # l 


2' 


l J 


2d 2 -5d(d+l) — -^<d+l) 
cr 


n 

7 


> 




n 


n 


> 






d+(d+l) 






then if 



D>0 (2.15) has i real root and two conjugate imaginary roots; 
D=0 (2.15) has 3 real roots, at least two of which are equal; 
D<0 (2.15) possesses 3 real and unequal roots. 



Note that if d=°° so that e n has a normal distribution, then certainly D>0 

U J 

and (2.2) has a unique maximum. If d<°° and C o n is small enough, then D>0 
and once again (2.2) will have a unique maximum. If d< oo and C*a n 2 is large 



enough (actually, larger than 
2 



(6V3-10)d 4d , 
d+1 ~d+L’ 



then D<0 for an interval of 



values of (y-m n .!) and (2.15) will have 3 real unequal roots; in this case (2.2) will 
have two local maxima. 



The likelihood procedure computes D. If D>0 it uses the ALMA 
procedure with weight 



to compute 0 n . If D<0, then two candidate estimates 0 b and 02 of 0 are 
computed. Both estimates are obtained via the ALMA procedure (2.7)-(2.9). 
One approximates weight (2.9) by setting as in (2.10); think of the result 

as prior-dominated. The other approximates weight (2.9) by setting 0=y,so that 



evaluated at each value of 0:0! and 0 2 The quoted estimate of 0 n is set equal to 
the 0, that comes closest to maximizing the global likelihood; the estimate of the 
variance is set equal to the corresponding C n . 

3. THE KALMAN AND WEST PROCEDURES. 

In this subsection, the traditional KALMAN procedure will be described 
for the model (1.1)-(1.2). A procedure proposed by West (1981) will also be 
discussed. 

3.1 The KALMAN Procedure. 

The KALMAN filter finds the estimate of 0 n which minimizes the 
conditional mean square error of (^ n *®n) given y n . If {e n } are independently 

normally distributed with mean 0 and variances {y n }» then the KALMAN filter 
can be viewed as a Bayesian updating procedure; see Meinhold and Singpurwalla 
(1983). 

The Bayesian KALMAN procedure assumes ©n-ily"' 1 is normal with mean 
m,,.! and variance Cn.i. Thus, from (1.1) ©ply"- 1 is normal with mean nip.! and 
variance C^=C n .! +x n . From (1.2) 



-l 




(2.17) 



w(y) = ; the result is data-determinated. The likelihood function is then 




n 



P(0 n €d0,Y n 6 dy | y"* 1 ) = K exp 



Yn 



d0dy 



(3.1) 



(3.2) 




Thus 0 n |y n has a normal distribution with mean 




n 



(3.3) 



and variance 




(3.4) 



The estimate of 0 n given y n is then 

0 n = m n 



(3.5) 



and an estimate of the variance of 0 n is C n . 

Comparing (3.3)-(3.4) with (2.10)-(2.12) indicates that, if y n is close to 
m n _i, then the ALMA procedure will closely resemble the KALMAN. In 

2 

particular, if y n = o n and d— the 2 estimators are identical. However, if y n is 

far from m n _i, then the ALMA procedure will tend to discount that observation, 

relying on its estimate of 0 n _i to strongly influence its estimate of 0 n . This 
behavior implies that the ALMA procedure will be less quickly responsive to 
changes in the values of 0 n than will be the KALMAN. This is the price paid for 
robustness to outlying measurement errors: KALMAN treats all changes in 
observations as representative of structural ( 0 n ) changes; ALMA is more 
tentative. Of course ALMA may be tuned towards KALMAN by increasing the 
d-value. 

3.2 The West Procedure. 



West proposes an estimation procedure for 0 n given y n in the case in 

normal, West's procedure reduces to the KALMAN filter. 

Once again, assume 9 n _i I y 0 ' 1 is normal with mean m n _i and variance 
C n _j so that 0 n ly nl is normal with mean m n .]and variance C n =C n _i+ x n 



which the density p r ^ is symmetnc about 0. In the special case in which 



Pie GdG, Y n edy|Yi=yi Y„.i=y n -i) 



= K exp 



Vm n l ) 2 -7 + ln p P (y- 6 ) 



d0<V 



(3.6) 



a K exp 



,2 1 



(e-m„-i y 



^(e-m n -i) — +1 In P e (y-m n .j) + g(y-m n .i)(e-m n -i) - G(y-m n .O ^ 

C n ' 



(3.7) 



where a Taylor expansion provides 

g(u) = £ ^(u) (3.8) 

G(u) = \ <u) - <3 9) 



Completing the square in (3.7) results in 
P[0 n ede,YnGdy|Yi=y 1 ,...,Y n .i=y n . 1 ] 



R 



exp <> 



+ G(y-m n .i) 



(0- m n-l)*g(y-m n -i)| — r + G(y-m n -i) 

c: 



l-!\ 



2 ) 



(3.10) 



Hence, P{0 n ed0 I Y 1 =y 1 ,...,Y n =y n } is approximated by a normal 
distribution having mean 



m n = m n .i + C n g(y n -m n .i) 



(3.11) 



and variance 



-l 



C n = 



+ G (yn-m n .i) 



(3.12) 



In the special case in which e n has a Student-t distribution with d degrees 
of freedom and scale parameter o n , 



p e (u)=c(d)— 

n a n 



1 

a n 



-(d+1) 



(3.13) 




(3.14) 



(3.15) 



Since G(y n -m n -i) is playing the role of a variance in (3.10), but may 
become embarrassingly negative for large u. West suggests that it be replaced by 
max(0,G(y n -m n _i)); this step has been taken in the simulations that illustrate the 
various procedures proposed here. West suggests another possibility in West et 
al. (1985). 

4. A SIMULATION EXPERIMENT. 

All simulations were carried out on an IBM 3033 AP computer at the 
Naval Postgraduate School. Random numbers were generated using the 
LLRANDOMII random number package; cf. Lewis and Uribe (1981). 

For each replication of the simulation the model of (1.1)-(1.2) is generated 
for n=0, 00. In the simulations reported below { oo n } are iid normal with 



mean zero and variance one. For each replication, estimates of 0 n given y n 
are computed using each of the procedures described above. The data collected 
are the estimation error £) n - 0 n for n=25, 50, 75, 100 and the estimate of 
variance C n , n=25, 50, 75, 100. The number of independent replications is 



1000. 



Tables 1 and 2 report results of the KALMAN and ALMA procedures for 
simulations in which (e n ) are iid normal with mean zero and variance one. The 
ALMA procedure actually uses the incorrect measurement error model that 
{e n } are iid Student-t with d=3 degrees of freedom and variance equal to one. 
Results for the ALMA procedure are shown for weights as in (2.10), for k=1.0 
and k=0.25. The procedure was iterated 0, 1, and 2 times. 

Table 1 shows statistics of § n -0 n for n=25, 50, 75, 100. As anticipated, the 
KALMAN procedure which uses the correct (normal) model exhibits the 
smallest variance of & n - 0 n . The ALMA procedure with k=0.25 and 0 iterations 
and the ALMA procedure with k=l and 1 iteration have the smallest variances 
for the ALMA procedures. 



Table 2 exhibits the estimates of the variance of 0 n , namely C n , for the 
ALMA procedure for n=25, 50, 75, 100. The KALMAN estimate of the 
variance is the constant 0.618 for all of these n. This constant is the limiting 

solution to equation (3.4) with x n = Y n =l; that is, with C= C n 

n- 



1 

(!+T + 1 

a simple quadratic with appropriate solution 




The variance of § n -0 n for the KALMAN procedure in Table 1 is close to 
the calculated 0.618. 

The mean values of C n for the ALMA procedure with k=0.25 and 0 
iterations and k=l with 1 iteration are about half that of the corresponding 
variances of § n -0 n m Table 1 . 



Tables 3-4 report results for a simulation in which {e n } are iid Student-t 
with 3 degrees of freedom and variance equal to 1. Table 3 reports statistics of 
the estimation error, & n -0 n , for tie KALMAN, ALMA, Biweight, Likelihood, 
and West procedures. As usual, the KALMAN procedure assumes {e n } are iid 

normal with mean 0 and variance 1. The other procedures assume {e n } are iid 
Student-t with 3 degrees of freedom and variance equal to 1. The ALMA 
procedure with k=.25 and no iterations exhibits the smallest variance of 9 n - 0 n . 
The more complicated Like’ihood procedure with k=0.25 and no iterations 
exhibits the next-smallest variance The ALMA with k=l and 1 iteration 
exhibits the third smallest variance. 

The Biweignt procedure implemented with the constants in the weight 
(2.13) a=5,7,9 and k-,23 and 1 , the procedure was iterated 0,1, and 2 times. 

The results for a=5 we e much were than those for a~7 and 9 indicating that 
a=5 is not large enough l suppress outlying values; they are not reported. 
Iterating the biweight p ? e jure 1 and 2 times did not improve the results for 
any values of a. The result- of I able 7 indicate that the bi weight procedure with 
the smallest variance re*, k-1 f' an i a--7 with no iterations. 

The West procedure deescribed ir» West (1981) as currently implemented 
does not do as w'ell as tine KALMAN. The statistics of C n in Table 4 seem to 
indicate that the difficulty 'S with the estimate of variance, C n ; the fix for 
negative G(y-m., j) makes it possible for C n to increase by one in successive 
times over long periods <. .' time. 

Table 4 exhibits die statistics of C n . The KALMAN procedure, the 
ALMA procedure with k~0.25 and 0 iterations, the ALMA procedure with k=l 
and 1 iteration, the Likelihood procedure with k=0.25 and 0 iterations and the 

Bi weight with k=l , a=7 all have mean C n approximately half the variance of £) n - 
9n- 

5. CONCLUSIONS. 



The simulation result; T anaed to date ndicate that a satisfactory robust 
KALMAN=ALM A f r cedun utilizes ±e k=0.25 weight-starting option and 
requires no iteration. W fine tde a >C\e filter is about 1% less efficient than the 
KALMAN when measurement errors are ideally Gaussian, it is about 6% more 
efficient when errors am long-tailed non-Gaussian; efficiency is in terms of 
ratios of (estimated) variances and is not the only meaningful criterion. 



Examination of Table 3 reveals through values of skewness, and kurtosis , that as 
anticipated, the robust ALMA estimation errors are substantially more closely 
Gaussian than are the corresponding KALMAN products when measurement 
errors are Student-t. 
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Table 1 

Statistics of § n -0 n 

Normal Measurement Errors with Variance 



Time n: 


25 


50 


75 


Proc Nbr k 


M V S K 


M V S K 


M V S K 


Iter 








K 


0.00 0.61 -0.04 0.01 


0.02 0.62 -0.01 -0.15 


-0.03 0.63 0.02 -0.29 


A 0 1.0 


-0.02 0.91 -0.07 0.60 


0.03 0.77 0.21 1.40 


-0.01 0.78 0.01 0.18 


0.25 


0.00 0.65 -0.02-0.07 


0.03 0.65 0.04 0.03 


-0.02 0.64 0.01 -0.14 


A 1 1.0 


0.01 0.70 -0.05 -0.02 


0.04 0.69 0.15 0.62 


-0.01 0.68 -0.00 0.02 


025 


0.02 0.70 0.00-0.07 


0.03 0.76 0.02 -0.09 


0.03 0.74 0.04-0.10 


A 2 1.0 


0.01 0.71 0.02-0.09 


0.04 0.75 0.06 0.06 


-0.02 0.72 0.01 -0.04 


0.25 


0.02 0.77 -0.01 -0.11 


0.02 0.83 0.02-0.10 


-0.03 0.82 0.04 -0.09 



Statistics 
M = Mean 
V = Variance 
S = Skewness 
K = Kurtosis 



1 



100 

M V S K 

0.00 0.60 -0.05 -0.19 
0.01 0.84 -0. 19 0.78 
0.00 0.67 -0.05 -0.17 
0.03 0.73 0.06 0.01 
0.00 0.75 0.02 030 
0.00 0.76 0.01 0.25 
0.00 0.81 0.04 031 



Procedure (Proc-i 
K = KALMAN 
A = ALMA 



Table 2 
Statistics of C n 

Normal Measurement Errors with Variance 1 



Time n: 


25 


50 


75 


100 


Proc Nbr 


k 


M 


V 


M 


V 


M 


V 


M 


V 


Iter 




















A 0 


1.0 


.50 


.08 


.49 


.07 


.48 


.06 


.48 


.07 




0.25 


.31 


.01 


.30 


.01 


.30 


.01 


.30 


.01 


1 


1.0 


.23 


.02 


.22 


.02 


.22 


.02 


.22 


.02 




0.25 


.04 


.00 


.14 


.00 


.14 


.00 


.13 


.00 


2 


1.0 


.14 


.01 


.13 


.01 


.13 


.01 


.13 


.01 




0.25 


.09 


.00 


.09 


.00 


.09 


.00 


.09 


.00 



Procedure (Proc.) 
A = ALMA 



Statistics 
M=Mean 
V = Variance 



Table 3 

Statistics of § n -0 n 

Student-t Measurement Errors with 3 degrees of freedom and Variance 1. 



Time n: 


25 


50 


75 


100 


Proc Nbr 


k a 


M 


V 


S K 


M V S K 


M V S K 


M V S K 




Iter 
















K 


- 


- 


0.01 


0.57 


- 0 . 4 $ 27 


0.02 0.53 0.02 22 


0.02 0.67 0.78 8.7 


0.04 0.54 - 0.17 1.7 


A 


0 


1.0 - 


0.03 


0.67 


- 0.07 1.0 


0.02 0.58 0.08 1.1 


0.02 0.71 0.02 1.4 


0.01 0 . 65 - 0.10 1.7 






0.25 - 


0.01 


0.53 


- 0.16 1.6 


0.01 0.48 0.08 \5 


0.02 037 0.01 21 


0.02 0.50 - 0 . 2 b 23 


A 


1 


1.0 - 


0.01 


0.55 


- 0.09 \5 


0.01 0.49 0.09 12 


0.02 0.61 0.01 1.8 


0.01 0 . 52 - 0.23 2.4 






0.25 - 


- 0.01 


0.63 


- 0.46 4.1 


0.01 0.58 0.03 3.1 


0.03 0.69 0.06 33 


0.03 0.58 0.20 2.7 


A 


2 


1.0 - 


- 0.01 


0.58 


- 0.16 2.4 


0.01 0.54 0.05 1.6 


0.03 0.64 0.02 1.9 


0.02 0.54 - 0.26 26 






0.25 - 


- 0.02 


0.71 


- 0.66 5_5 


0.01 0.66 0.09 43 


0.03 0.79 0.18 4.9 


0.04 0.64 - 0.05 29 


B 


0 


1.0 7 


0.01 


0.56 


- 0.17 25 


0.01 031 0.07 22 


0.02 039 0.09 2.6 


0.02 0 . 52-031 24 






025 7 


0.01 


0.61 


- 0.68 52 


0.01 0.57 0.04 3.9 


0.03 0.69 0.35 6.0 


0.04 0 . 56 - 0.13 2.7 


B 


0 


1.0 9 


- 0.01 


0.57 


- 0.42 3.5 


0.01 034 0.00 3.0 


0.02 0.62 0.07 3.3 


0.03 0.53 - 0.22 23 






0.25 9 


- 0.01 


0.63 


• 0.76 5.7 


0.01 037 0.05 4.1 


0.03 0.72 0.75 9.9 


0.04 0 . 57 - 0.11 28 


L 


0 


1.0 - 


0.03 


0.68 


0.12 \5 


0.03 036 0.09 1.0 


0.03 0.69 0.02 13 


0.01 0 . 65 - 0.13 1.8 






025 - 


0.01 


054 


0.23 1.7 


0.01 0.48 0.11 13 


0.03 0.60 0.04 21 


0.02 0.52 0.29 25 


L 


1 


1.0 - 


0.01 


055 


0.13 1.7 


0.01 0.49 0.09 12 


0.02 0.61 0.01 1.8 


0.01 0.52 0.23 24 






025 - 


- 0.01 


0.63 


0.53 4.1 


0.00 037 0.07 28 


0.03 0.69 0.07 3.4 


0.03 0.58 0.20 27 


L 


2 


1.0 - 


0.00 


0.59 


0.22 25 


0.01 033 0.08 1.6 


0.03 033 0.02 1.9 


0.02 0.54 - 0.26 26 






025 - 


- 0.02 


0.72 


0.70 5.5 


0.00 0.65 0.01 3.9 


0.03 0.79 0.18 4.9 


0.04 0.64 - 0.05 29 


W 


- 


- - 


- 0.06 


112 


026 3.2 


031 384 0.14 43 


0.08 774 0.03 5.1 


- 0 . 791249 - 0.25 5.2 



Procedure (Troc) 
K = KALMAN 
A = ALMA 
B = Biweight 
L = Likelihood 
W = West 



Statistics 
M = Mean 
V = Variance 
S = Skewness 
K = Kurtosis 



Table 4 
Statistics of C n 

Student-t Measurement Errors with 3 degrees of freedom and Variance 1. 



Time n: 


25 


50 


75 


100 


Proc 


Nbr k 
Iter 


a 


M 


V 


M 


V 


M 


V 


M 


V 


A 


0 1.0 




.46 


.07 


.44 


.06 


.48 


.07 


.45 


.06 




0.25 


- 


.29 


.01 


.28 


.01 


.29 


.01 


.29 


.01 


A 


1 1.0 


- 


.21 


.02 


.20 


.02 


.21 


.02 


.21 


.02 




0.25 


- 


.13 


.00 


.13 


.00 


.14 


.00 


.14 


.00 


A 


2 1.0 


- 


.13 


.01 


.12 


.01 


.13 


.01 


.13 


.01 




0.25 


. 


.09 


.00 


.09 


.00 


.09 


.00 


.09 


.00 


B 


0 1.0 


7 


.29 


.01 


.29 


.00 


.29 


.01 


.29 


.00 




0.25 


7 


.27 


.00 


.27 


.00 


.27 


.00 


.27 


.00 


B 


0 1.0 


9 


.28 


.00 


.28 


.00 


.28 


.00 


.28 


.00 




0.25 


9 


.27 


.00 


.27 


.00 


.27 


.00 


.27 


.00 


L 


0 1.0 


- 


.44 


.06 


.44 


.06 


.46 


.06 


.46 


.06 




0.25 


- 


.29 


.01 


.28 


.01 


.29 


.01 


.29 


.01 


L 


1 1.0 


- 


.21 


.02 


.20 


.02 


.21 


.02 


.21 


.02 




0.25 


- 


.14 


.00 


.14 


.00 


.14 


.00 


.14 


.00 


L 


2 1.0 


- 


.13 


.01 


.12 


.01 


.13 


.01 


.13 


.01 




0.25 


- 


.09 


.00 


.09 


.00 


.09 


.00 


.09 


.00 


W 


- - 


- 


8.8 


64 


16 


240 


23 


543 


29 


946 



Procedures (Pr oc. 'i 
A = ALMA 
B = Biweight 
L= Likelihood 
W = West 



Statistics 
M=Msan 
V= Variance 
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